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Abstract 



Universality of eigenvalue spacings is one of the basic characteristics of random 
matrices. We give the precise meaning of universality and discuss the standard 
CN ! universality classes (sine, Airy, Bessel) and their appearance in unitary, orthog- 

| onal, and symplectic ensembles. The Riemann-Hilbert problem for orthogonal 

£S) ' polynomials is one possible tool to derive universality in unitary random matrix 

ensembles. An overview is presented of the Deift/Zhou steepest descent analy- 
sis of the Riemann-Hilbert problem in the one-cut regular case. Non-standard 
universality classes that arise at singular points in the spectrum are discussed 
at the end. 



6.1 Heuristic meaning of universality 



Universality of eigenvalue spacings is one of the basic characteristic features 
of random matrix theory. It was the main motivation of Wigner to introduce 
random matrix ensembles to model energy levels in quantum systems. 

On a local scale the eigenvalues of random matrices show a repulsion which 
implies that it is highly unlikely that eigenvalues are very close to each other. 
Also very large gaps between neighboring eigenvalues are unlikely. This is in 
sharp contrast to points that are sampled independently from a given distribu- 
tion. Such points exhibit Poisson statistics in contrast to what is now called 
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GUE/GOE/GSE statistics. These new random matrix theory type statistics 
have been observed (and sometimes proved) in many other mathematical and 
physical systems as well, see |Dei07d] . and are therefore universal in a sense 
that goes beyond random matrix theory. Examples outside of random matrix 
theory are discussed in many chapters of this handbook. In this chapter the 
discussion of the universality is however restricted to random matrix theory. 

The universality conjecture in random matrix theory says that the local 
eigenvalue statistics of many random matrix ensembles are the same, that is, 
they do not depend on the exact probability distribution that is put on a set 
of matrices, but only on some general characteristics of the ensemble. The 
universality of local eigenvalue statistics is a phenomenon that takes place for 
large random matrices. In a proper mathematical formulation it is a statement 
about a certain limiting behavior as the size of the matrices tends to infinity. 

The characteristics that play a role in the determination of the universality 
classes are the following. 

• Invariance properties of the ensemble, of which the prototype is invariance 
with respect to orthogonal, unitary, or unitary-symplectic conjugation. 
As discussed in Chapters 4 and 5 of this handbook, these lead to random 
matrix ensembles with an explicit joint eigenvalue density of the form 

1 n 
p(xi,...,ar n ) = — ]"[ \xj-x^T\e- v ^ (6.1.1) 

" l<i<j<n j=l 

where /3 = 1,2,4 corresponds to orthogonal, unitary, and symplectic en- 
sembles, respectively. The case V(x) = \x 1 gives the Gaussian ensembles 
GOE, GUE, and GSE. The local eigenvalue repulsion increases with j3 
and different values of (3 give rise to different universality results. 

• An even more basic characteristic is the distinction between random ma- 
trix ensembles of real symmetric or complex Hermitian matrices. Even 
without the invariance assumption, the universality of local eigenvalue 
statistics is conjectured to hold. This has been proved recently for large 
classes of Wigner ensembles, i.e., random matrix ensembles with inde- 
pendent, identically distributed entries, see |ErdlOa[ lErdlOb} ITaoll] and 
the survey [ErdlOcj . See Chapter 21 of this handbook for more details. 
Universality is also expected to exist in classes of non- Hermitian matrices, 
see Chapter 18. 

• Another main characteristic is the nature of the point around which the 
local eigenvalue statistics are considered. A typical point in the spectrum 
is such that, possibly after appropriate scaling, the limiting mean eigen- 
value density is positive. Such a point is in the bulk of the spectrum, and 
universality around such a point is referred to as bulk universality. 
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At edge points of the limiting spectrum the limiting mean eigenvalue 
density typically vanishes like a square root, and at such points a different 
type of universality is expected, which is known as (soft) edge universality. 
At natural edges of the spectrum, such as the point zero for ensembles of 
positive definite matrices, the limiting mean eigenvalue density typically 
blows up like the inverse of a square root, and universality at such a point 
is known as hard edge universality. 

• At very special points in the limiting spectrum the limiting mean eigen- 
value density may exhibit singular behavior. For example, the density 
may vanish at an isolated point in the interior of the limiting spectrum, 
or it may vanish to higher order than square root at soft edge points. This 
non-generic behavior may take place in ensembles of the form (|6.1.ip with 
a varying potential NV 



as n, N — y co with n/N — > 1. 

At such special points one expects other universality classes determined 
by the nature of the vanishing of the limiting mean eigenvalue density at 
that point, see also a discussion in Chapter 14. 

In the rest of this chapter we first give the precise meaning of the notion 
of universality and we discuss the limiting kernels (sine, Airy and Bessel) as- 
sociated with the bulk universality and the edge universality for the unitary, 
orthogonal, and symplectic universality classes. In Section [6.31 we discuss uni- 
tary matrix ensembles in more detail and we show that universality in these 
ensembles comes down to the convergence of the properly scaled eigenvalue 
correlation kernels. In Section 16.41 we discuss the Riemann-Hilbert method in 
some detail. The Riemann-Hilbert method is one of the main methods to prove 
universality in unitary ensembles. In the final Section 16.51 we discuss certain 
non-standard universality classes that arise at singular points in the limiting 
spectrum. We describe the limiting kernels for each of the three types of sin- 
gular points, namely interior singular points, singular edge points, and exterior 
singular points. 

Other approaches to universality are detailed in Chapter 21 for Wigner 
ensembles and in chapter 16 using loop equation techniques. 

6.2 Precise statement of universality 

For a probability density p n (x±, . . . ,x n ) on R n that is symmetric (i.e., invari- 
ant under permutations of coordinates), let us define the /c-point correlation 




n 
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function by 

Rn t k(xi,. ..,x k )= ^— - J ■■■ J p n (x 1 ,...,x n )dx k +i---dx n . (6.2.1) 

Up to the factor n!/(n — k)\ this is the kth marginal distribution of p n . 

Fix a reference point x* and a constant c n > 0. We center the points around 
x* and scale by a factor c n , so that (xi, . . . , x n ~) is mapped to 

(c n (x\ x ),..., c n (x n X )). 

These centered and scaled points have the following rescaled /c-point correlation 
functions 

±Rn, k (x* + x* + ^, . . . , x* + ^ V (6.2.2) 

The universality is a property of a sequence (p n ) of symmetric probability 
density functions. Universality at x* means that for a suitably chosen sequence 
(c n ) the rescaled /c-point correlation functions (|6.2.2p have a specific limit as 
n — > oo. The precise limit determines the universality class. 

As discussed in chapter 4 for determinantal point processes, the main spec- 
tral statistics such as gap probabilities and eigenvalue spacings can be expressed 
in terms of the /c-point correlation functions. If the limits (|6.2.2|) exist and be- 
long to a certain universality class, then this also leads to the universal behavior 
of the local eigenvalue statistics. 

6.2.1 Unitary universality classes 

The unitary universality classes are characterized by the fact that the limit of 
(|6.2.2p can be expressed k determinant 

det [Kixuxj)]^^ 

involving a kernel K(x,y), which is called the eigenvalue correlation kernel. 

Bulk universality A sequence of symmetric probability density functions 
(p n ) then has bulk universality at a point x* if there exists a sequence (c n ), so 
that for every k, the limit of (|6.2,2p is given by 

det ; /r n, i.r, .-•,);, , r , 

where 

K^x,y) = S ^^l (6.2.3) 
tt{x - y) 

is the sine kernel. Bulk universality is also known as GUE statistics. 
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Edge universality The (soft) edge universality holds at x* if the limit of 
(16X21) is equal to 

det [ir Ai (^,x i )] 1 < iJ < fe 

where 

rM = A\(x) A\'{y) - Ai'(x) Ai(y) 

x-y 

is the Airy kernel. Here Ai denotes the Airy function, which is the solution of 
the Airy differential equation y"(x) = xy(x) with asymptotics 

Ai(x) = — — , lA e~^ x3/2 (l + 0{x~ 3/2 )) , as x -> +oo. 
2 v / 7rx 1 / 4 V / 

The Airy kernel is intimately related to the Tracy- Widom distribution for the 
largest eigenvalue in random matrix theory |Tra94j . 



Hard edge universality A hard edge is a boundary for the eigenvalues that 
is part of the model. For example, if the random matrices are real symmetric 
(or complex Hermitian) positive definite then all eigenvalues are non-negative 
and zero is a hard edge. 

The hard edge universality holds at x* if the limit of (|6.2.2p is equal to 

det [K^^x,)]^^ 

where 

2(x - y) 

and J a is the usual Bessel function of order a. 

The Bessel kernels depend on the parameter a > — 1 which may be inter- 
preted as a measure of the interaction with the hard edge. The bigger a, the 
more repulsion from the hard edge. 



6.2.2 Orthogonal and symplectic universality classes 

The orthogonal and symplectic universality classes are characterized by the fact 
that the limit of (|6.2.2|) is expressed Pfaffian 

PflKfaxj)]^^ (6.2.6) 

of a 2 x 2 matrix kernel 

K(x,y)=(^ y \ f 12 S X ' y |V (6.2.7) 
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Recall that the Pfaffian Pf(A) of a 2k x 2k skew symmetric matrix A is such 
that 

Pi(A) = y/det(A). 

The matrix in ()6.2.6p is a 2k x 2k matrix written as a x fc block matrix with 
2x2 blocks. It is skew symmetric provided that K(y, x) = —K(x, y) T where T 
denotes the matrix transpose. 



Bulk universality A sequence of symmetric probability density functions 
(p n ) has orthogonal/symplectic bulk universality at a point x* if there exists a 
sequence (c n ), so that for every k, the limit of ()6.2.2|) is given by 



Pf 



K '^(x^, Xj ) 



l<i,j<k 



with 



K 



sin,fl=l 



/ d sin7r(x — y) 
dx ir(x — y) 
sin7r(x — y) 

V 



sin7r(x — y) 
ir(x - y) 



x y sin 7rf 



-dt 



1 



•sg: 



n(x - y) 



vr(x - y) J 7rt 2 

where sgn(t) = 1, 0, —1 depending on whether t > 0, t = 0, or t < 0, in the case 
of orthogonal (i.e., (5 = 1) bulk universality, and 

/ d sin2-7r(x — y) sin2-7r(x — y) \ 
dx 2ir(x — y) 2tt(x — y) 

sin27r(x — y) f x ~ y sin27rt 
Jo 



K Bin ' P= \x,y) 



\ 



-dt 



2ir(x-y) J 2nt ) 
in the case of symplectic (i.e., /3 = 4) bulk universality. 



Edge universality The orthogonal/symplectic (soft) edge universality holds 
at x* if the limit of ()6.2.2j) is equal to 



Pf 



Xi , Xj j 



l<i,j<k 



with 



K 



Ai,/3=1 
11 



^-K K \x,y) + l -k\{x)A\{y) 



K^= l (x,y) = K Al (x,y) + i ; Ai(x)- 



Ai(t) dt 



K. 



Ai,/3=1 
21 1 

Ai,/3=1 
22 



x,y) 



-K 



x,y) 



Ai,/3=1 
12 1 

K A[ (t,y)dt-- / M(t)dt 



roc roo | 

+ -J Ai(t) dt ■ J Ai(t) dt - - sgn 



(6.2.8) 



(x - y) 
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in the case of orthogonal edge universality, and 

K^=\x,y) = ±JLK M (x,y) + \ Ai(x) Ai(y) 

K^=\x, y) = -K Ai (x, y)-- Ai(x) • / Ai(i) di 

2 4 A (6.2.9) 

K A ^ =4 (x,y) = -K A ^ =4 (x,y) 

-| /"OO I /"OO /'OO 

K^ = \x,y) = --J K Ai (t,y)dt + -j A\{t)dt.j M(t)dt 

in the case of symplectic edge universality. The kernel if Al that appears in 
flgXgp and (|6T9"j) is the Airy kernel from ([5^|> . 

Hard edge universality The orthogonal/symplectic hard edge universality 
is expressed in terms of 2 x 2 kernels with Bessel functions. See Forrester [For 10] 
for the precise statement. 

6.2.3 Determinantal and PfafRan point processes 

The universality limits have the characteristic properties of determinantal or 
Pfafnan point processes, see [ForlOl ISosOOj and Chapter 11 of this handbook. 

If the probability densities p n on M n also arise from determinantal point pro- 
cesses, then the statement of universality can be expressed as the convergence 
of the corresponding correlation kernels after appropriate scaling. 

The probability density p n arises from a determinantal point process, if 
there exist (scalar) kernels K n so that for every n and k, 

Rn,k{xii ■ ■ ■ > Xfi) = det \K n \Xii ■Ej)]l<i j<k • 
In particular, one then has for the 1-point function (particle density), 

R n ,i(x) = K n (x,x) 
and for the probability density p n itself 

p n (xi,...,x n ) = det [Knix^Xj)}^^. 

In this setting the bulk universality comes down to the statement that the 
centered and rescaled kernels 

—K n fx* + —,x* + — ^ (6.2.10) 



tend to the sine kernel (|6.2.3j) as n — > oo, and likewise for the edge universality. 
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Similarly, the probability densities p n come from Pfaffian point processes, 
if there exist 2x2 matrix kernels K n so that for every n and k, 

Rn,k ; ■ ■ ■ > % k) — Pf \K n (xj , Xj )] ,j<k ' 

A proof of universality in a sequence of Pfaman point processes then comes 
down to the proof of a scaling limit for the matrix kernels as n — > oo. 

Unitary random matrix ensembles are basic examples of determinantal point 
processes, while orthogonal and symplectic matrix ensembles are examples of 
Pfaman point processes. 

For the classical ensembles that are associated with Hermite, Laguerre and 
Jacobi polynomials the existence and the form of the limiting kernels has been 
proven using the explicit formulas that are available for these classical orthog- 
onal polynomials, see e.g. [For 101 lMeh04] . 

For non-classical ensembles, the results about universality are fairly com- 
plete for unitary ensembles, due to their connection with orthogonal polyno- 
mials. This will be discussed in more detail in the next section, and we give 
references there. 

The first rigorous results on bulk universality for orthogonal and symplectic 
ensembles are due to Stojanovic [StoOO. who discussed ensembles with a quartic 
potential. This was extended by Deift and Gioev [Dci07a] for ensembles (]6.1.ip 
with polynomial V. Their techniques are extended to treat edge universality in 
[Dei07bl IDei07c| , see also the recent monograph [De i09] . Varying weights are 
treated by Shcherbina in |Shc08[ IShc 09a] . 

6.3 Unitary random matrix ensembles 

We explain in more detail how the universality classes arise for the eigenvalues 
of a unitary invariant ensemble 

J_ e ~NTrV(M) dM (631) 
Z n ,N 

defined on the space of n x n Hermitian matrices, where 

n 

dM = Y[ dMa ] [ d Re d Im My , 

1=1 l<i<j<n 

and Z n< N is a normalization constant, see also Chapter 4. The potential function 
V in ()6.3.ip is typically a polynomial, but could be more general as well. To 
ensure that (|6.3.ip is well-defined as a probability measure, we assume that 



(6.3.2) 
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The factor N in (|6.3.ip will typically be proportional to n. This is needed in 
the large n limit in order to balance the repulsion among eigenvalues due to the 
Vandermonde factor in f|6.3.3[) and the confinement of eigenvalues due to the 
potential V. 

6.3.1 Orthogonal polynomial kernel 

The joint probability density for the eigenvalues of a matrix M of the ensemble 
()6.3.ip has the form 

1 " 

Z n M 

' l<i<j<n j=l 

Introduce the monic polynomials jv, Pk,N{x) = x k + ■ ■ ■ , that are orthogonal 
with respect to the weight e~ NV ( x ) on R: 



0, for j = 0,...,k-l, 

and the orthogonal polynomial kernel 



P k)N {x)x^e- NV ^dx=r 2 n L ; 1J _ X ' (6.3.4) 



K n , N (x, y) = V^V^) Pk ' N{x) 2 Pk ' Niy) . (6.3.5) 

k=0 1 k,N 

Using the formula for a Vandermonde determinant and performing elementary 
row operations we write (|6.3.3p as 



^ 2 

Pu,n{xi, ... ,x n ) = (det [Pk-iM x j)]i<j,k<n) II 



n 

e -NV(x 3 ) 

3=1 



l\k=0 7fc,AT 



Z 



det 



% / c ~NV(xA P k~l,N(Xj) 
lk-l,N 



l<j,k<n / 



Evaluating the square of the determinant using the rule (det A) 2 = det(AA T ) 
we obtain that (|6.3.3[) can be written as 

p n ^ N (x 1 , . . . ,x n ) = — detlK^N^XijXj)]^^^, (6.3.6) 

since it can be shown that Z U) n = n\ IIfc=o 7fc at- The orthogonality condition 
(|6,3.4p is then used to prove that the fc-point correlation functions (|6.2.ip are 
also determinants 

Rn,k{xi, ... ,x k ) = det [K n ^ N (x i ,Xj)] 1<i j<k , 
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which shows that the eigenvalues are a determinantal point process with corre- 
lation kernel K n jsj. For the above calculation, see also Chapter 4. 

The eigenvalues of orthogonal and symplectic ensembles of random matrices 
have joint probability density 



n,JV l<i<j<n 3=1 



-NV(xj) 



with /3 = 1 for orthogonal ensembles and (3 = 4 for symplectic ensembles. These 
probability densities are basic examples of Pfaffian ensembles, as follows from 
the calculations in e.g. [DeiO.91 [ Tra98j , see also Chapter 5 of this handbook. 

6.3.2 Global eigenvalue regime 

As n, N — > oo such that n/N — > 1 the eigenvalues have a limiting distribution. 
A weak form of this statement is expressed by the fact that 

lim — K n ^(x,x) = pv(x), i£l (6.3.7) 
n,AT->oo n ' 
n/N-+l 

exists. The density py describes the global or macroscopic eigenvalue regime. 
The probability measure dpy(x) = py{x)dx minimizes the weighted energy 

log , . du(x)du(y) + I V(x)du(x) (6.3.8) 
\x-y\ J 

among all Borel probability measures p on R. Heuristically, it is easy to un- 
derstand why (|6.3.8|) is relevant. Indeed, from (|6.3.3|) we find after taking log- 
arithms and dividing by — n 2 , that the most likely distribution of eigenvalues 
xi, . . . ,x n for (|6.3.3p is the one that minimizes 

i,j=i j=i 

This discrete minimization problem leads to the minimization of (|6.3.8j) in the 
continuum limit as n, N — > oo with n/N -4 1. 

The minimizer of (|6.3.8|) is unique and has compact support. It is called 
the equilibrium measure in the presence of the external field V, because of 
its connections with logarithmic potential theory [Saf97| . The proof of (|6.3.7p 
with py(x)dx minimizing (|6.3.8p is in [Dei99a, Joh98j. See also the remark 
at the end of subsection 16.4.81 below. It follows as well from the more general 
large deviation principle that is associated with the weighted energy (|6.3.8p , see 
[Ben97| and also [AndlOl §2.6]. See also Chapter 14 of this handbook. 
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If V is real analytic on R then py is supported on a finite union of intervals, 

say 

m 

suppOy ) = |J [aj,bj], 
i=i 

and the density py takes the form 

pv(x) = —hj(x)\l (bj — x)(x — cij), x € [dj, bj], j = l,...,m (6.3.9) 

7T V 

where hj is real analytic and non-negative on [a,-, bj], see |Dei98] . Another 
useful representation is that 

py{x) = ^q v (x), xeR, (6.3.10) 

where q v is the negative part of 

If V is a polynomial then qy is a polynomial of degree 2(degV — 1). In that 
case the number of intervals in the support is bounded by ^degV, see again 
[Dei98]. See also a discussion in Chapters 14 and 16. 

6.3.3 Local eigenvalue regime 

In the context of Hermitian matrix models, the universality may be stated as the 
fact that the global eigenvalue regime determines the local eigenvalue regime. 
The universality results take on a different form depending on the nature of the 
reference point x*. 



Bulk universality A regular point in the bulk is an interior point x* of 
supp( / uy) such that py{x*) > 0. At a regular point in the bulk one has 

lim — K n , N (x* + — , x* + 1L) = K s ' m (x, y) (6.3.12) 
n^oo cn \ cn cn) 

where c = py{x*) and K sm is the sine kernel (|6.2.3p . 

Convincing heuristic arguments for the universality of the sine kernel were 
given by Brezin and Zee [Bre93 by the method of orthogonal polynomials. 
Supersymmetry arguments were used in |Hac95] . Rigorous results for bulk 
universality (|6.3.12p beyond the classical ensembles were first given by Pastur 
and Shcherbina |Pas97| . and later by Bleher and Its [Ble99] and by Deift et 
al. |Dei99bl IDei99c| for real analytic V. In these papers the Riemann-Hilbert 
techniques were introduced in the study of orthogonal polynomials, which we 
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will review in section [6T4"l below Conditions on V in (|6,1.2p were further weakened 
in |McL08llPas08] . 

In more recent work of Lubinsky [Lub08a, Lub09a, Lub09b] and Levin and 
Lubinsky [Lev08, Lev09] it was shown that bulk universality holds under ex- 
tremely weak conditions. One of the results is that in an ensemble (|6,l.ip 
restricted to a compact interval, bulk universality holds at each interior point 
where V is continuous. 



Soft edge universality An edge point x* 6 {a\, bx, ... , a m , b m } is a regular 
edge point of supply) if hj(x*) > in (I6.3.9j) . In that case the density py 
vanishes as a square root at x* . The scaling limit is then the Airy kernel (|6.2.4I) . 
If x* = bj is a right-edge point then for a certain constant c > 0, 



n— s-oo 



while if x* = aj is a left-edge point, we find the same limit after a change of 
sign 



^ r^m K "," x *-r^73' x *-Hw3 = * 0^)- ( 6 - 3 - 14 ) 

Thus the regular edge points belong to the Airy universality class. 

The Airy kernel was implicitly derived in |Bow91j . For classical ensembles 
the soft edge universality (|6.3.13p is derived in [For93l Nag93| , see also |ForlO| . 



For quartic and sextic potentials V, it is derived in [Kan97| by the so-called 
Shohat method. For real analytic potentials potentials it is implicit in |Dei99c| 
and made explicit in [De i07b| . 



Hard edge universality The Bessel universality classes arise for eigenvalues 
of positive definite matrices. Let a > —1 be a parameter and consider 

J— (det M) a e- NT " v ^dM 

Z n ,N 

as a probability measure on the space of n x n Hermitian positive definite 
matrices. For the case V(x) = x this is the Laguerre Unitary Ensemble, which 
is also known as a complex Wishart ensemble, see Chapter 4. 

Assuming that a remains fixed, the global eigenvalue regime does not de- 
pend on a. The eigenvalue density py lives on [0, oo), and for x G [0, oo) it 
continues to have the representation (I6.3.10P but now qy is modified to 

Qv(x) = (pf) 2 ~ J V ' iX) x ~_ V s ' {S) dMs) ~\j V'(s)dpy(s). (6.3.15) 
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If J V'(s)dfjLv(s) > then is in the support of p,y and py has a square-root 
singularity at x = 0. 

The effect of the parameter a is noticeable in the local eigenvalue regime 
near x* = 0. If J V '(s)d//y(s) > 0, then the limiting kernel is the Bessel kernel 
(I6.2.5P of order a. For x, y > and for an appropriate constant c > 0, we have 

hm -L^ K n>N (j^, - JQ = ^ Bes ' a (x, y). (6.3.16) 

The hard edge universality (16.3, 16|) was proven in [Kui02, Lub08b]. 

Spectral singularity A related class of Bessel kernels 

- R ,- «/ Q+ l(vr2;)J Q _i (vry) - J i(vrx)J a+ i (vry) 
K Bes '<*(x, y) = n^Vv— 2 2{x _y ) 2 ^ (6.3.17) 

appears as scaling limits in Hermitian matrix models of the form 

J—\detM\ 2a e- NTlV ^dM (6.3.18) 

Z n ,N 

with a > -1/2. 

The extra factor |detM| 2a is referred to as a spectral singularity and the 
matrix model f|6.3. 18j) is relevant in quantum chromo dynamics, where it is re- 
ferred to as a chiral ensemble, see Chapter 32 of this handbook. The spectral 
singularity does not change the global density py of eigenvalues, but it does 
have an influence on the local eigenvalue statistics at the origin. Assuming that 
c = pv(0) > 0, one now finds 

hm — K n , N (—, JL) = K BeB > a (x, y), x,y>0, (6.3.19) 



n->oo cn \cn ctij 
instead of (|6.3.12|) . 

The universality of the Bessel kernels (I6.2.5P and (16.3. \1\ was discussed by 
Akemann et al. [|Ake_97b] for integer values of a. It was extended to non-integer 
a in [Kan98 . The model (16.3. 18D was analyzed in |Kui03| with the Riemann- 
Hilbert method. 

See }Kle00| for the spectral universality in orthogonal ensembles. 

Weight with jump discontinuity More recently |Foull| . a new class of 
limiting kernels was identified for unitary ensembles of the form 

^ e -TrV(M) dM 



Z, 



ii 



defined on Hermitian matrices with eigenvalues in [—1, 1] in cases where the 
weight function w(x) = e~ v ^ x \ x G [—1, 1], has a jump discontinuity at x = 0. 
The limiting kernels in [Foull| are constructed out of confluent hyper geometric 
functions, see also |Its08| . 
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6.4 Riemann-Hilbert method 

A variety of methods have been developed to prove the above universality results 
in varying degrees of generality. One of these methods will be described here, 
namely the steepest descent analysis for the Riemann-Hilbert problem (RH 
problem). This is a powerful method to obtain strong and uniform asymptotics 
for orthogonal polynomials. One of the outcomes is the limiting behavior of 
the eigenvalue correlation kernels. However, the method gives much more. It is 
also able to give asymptotics of the recurrence coefficients, Hankel determinants 
and other notions associated with orthogonal polynomials. In this section we 
closely follow the paper [Dei99c] of Deift et al, see also [Dei99a| . 

6.4.1 Statement of the RH problem 

A Riemann-Hilbert problem is a jump problem for a piecewise analytic function 
in the complex plane. The RH problem for orthogonal polynomials asks for a 
2x2 matrix valued function Y satisfying 

(1) Y : C \ M -> C 2x2 is analytic, 

(2) Y has boundary values on the real line, denoted by Y±(x), where Y + (x) 
(Y_(x)) denotes the limit of Y{z) as z — > x G R with Imz > (Imz < 0), 



The unique solution is given in terms of the orthogonal polynomials (|6.3.4p by 



The RH problem for orthogonal polynomials and its solution ([6.4. ip are due to 
Fokas, Its, and Kitaev |Fok92] . 

By the Christoffel-Darboux formula for orthogonal polynomials, we have 
that the orthogonal polynomial kernel ([6.3.5P is equal to 

K n ,N(x,y) = 



y^)^)^. N P nM*)Pn-l,N(y) ~ Pn-lM*)Pn,N(y) ( g A2) 



and 






(6.4.1) 



x-y 
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which in view of (|6.4.ip and the fact that det Y(z) = 1 can be rewritten as 

K n)N (x,y) = 

27Tii l_y/ e ~ NV{x) ^ e ~ NV{y) (° l)Y+(y)- l Y + (x)(^J, (6.4.3) 

for x, y € M. Other notions related to the orthogonal polynomials are also con- 
tained in the solution of the Riemann-Hilbert problem. The monic polynomials 
satisfy a three-term recurrence relation 

xP n>N (x) = P n+ i )N (x) + b n , N P n , N {x) + a n:N P n _ 1)N (x) 

with recurrence coefficients a n ^ > and b n ^N G R. The recurrence coefficients 
can be found from the solution of the RH problem by expanding Y around oo: 

as z — > oo. The 2x2 matrices Y\ and I2 do not depend on z, but do depend 
on n and N . Then 

a n ,N = (li)i a (^1)21 , ^,^v = - (ii) 22 ■ (6.4.4) 

( y l)l2 

6.4.2 Outline of the steepest descent analysis 

The steepest descent analysis of RH problems is due to Deift and Zhou |Dei93j . 
It produces a number of explicit and invertible transformations leading to a 
RH problem for a new matrix-valued function R with identitiy asymptotics at 
infinity. Also R depends on n and N, and as n,N — > 00 with n/N — > 1, the 
jump matrices tend to the identity matrix, in regular cases typically at a rate 
of 0(l/n). Then it can be shown that R(z) tends to the identity matrix as 
n, N — > 00 with n/N — > 1 and also at a rate of 0(l/n) in regular cases. 

Following the transformations in the steepest descent analysis we can then 
find asymptotic formulas for Y and in particular for the orthogonal polynomial 
P n ,N, the recurrence coefficients and the correlation kernel. 

With more work one may be able to obtain more precise asymptotic infor- 
mation on R. For example, if n = N and if we are in the one-cut regular case, 
then there is an asymptotic expansion for R(z): 



D , , T rW(z) R^{ 
R(z) ~ I + ^ + £ 



n n z 



with explicitly computable matrices R^\z). This in turn leads to asymptotic 
expansions for the orthogonal polynomials as well. We will not go into this 
aspect here. 
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Here we present the typical steps in the Deift-Zhou steepest descent analysis. 
We focus on the one-cut case, that is on the situation where V is real analytic 
and the equilibrium measure fj,y is supported on one interval [a, b]. We also 
assume that we are in a regular case, which means that 

pv{x) = ^Y^L = Lmx)W (b — x)(x — a), x G [a, b], 
ax ir 

with h{x) > for x G [a, b], and strict inequality holds in the variational 
condition (I6.4.6P below. Generically these regularity conditions are satisfied, 
see [KuiOO] . The singular cases lead to different universality classes, and this 
will be commented on in section 16.51 
For convenience we also take N = n. 



6.4.3 First transformation: normalization at infinity 

The equilibrium measure \iy is used in the first transformation of the RH prob- 
lem. The equilibrium measure satisfies for some constant I 

2 f log -J— + V{x) = £, xe[a,b], (6.4.5) 
J \x-y\ 

2 / log— —dfj, v (y) + V(x) > £, xeR\[a,b}. (6.4.6) 
J \x - y\ 

These are the Euler-Lagrange variational conditions associated with the min- 
imization of the weighted energy (|6.3.8[) . Since (|6.3.8p is a strictly convex 
functional on finite energy probability measures, the conditions (|6.4.5p -( |6.4.6p 
characterize the minimizer \iy . 

The equilibrium measure [iy leads to the ^-function 

g(z) = J log(z — x) dfiy(x), z G C \ (— oo, b], (6.4.7) 

which is used in the first transformation Y t— >• T. We put 

T(*)={ e ne/2)Y(z)[ e nG,W-</a)J- ( 6 - 4 - 8 ) 

The jumps in the RH problem for T are conveniently stated in terms of the 
functions 



( z )= I h(s){(s - b){s - a)) 1 ' 2 ds (6.4.9) 
h(s)((s-b)(s-a)) 1/2 ds, (6.4.10) 



where we use the fact that V is real analytic and therefore h has an analytic 
continuation to a neighborhood of the real line. Then T satisfies the RH problem 
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Figure 6.1: Contour Eg and jump matrices for the RH problem for S. This 
figure is reproduced from |Kui09| . 



(1) T is analytic in C \ R. 

(2) On R we have the jump T + = T_ Jt where 



( ' e 2n<t> + {x) 



Jt(x) 




(6.4.11) 



x < a. 



(3) T(z) = 1 + 0(l/z) as z -> oo. 

The RH problem is now normalized at infinity. The jump matrices on 
(— oo,a) and (b, oo) tend to the identity matrix as n — > oo, since <f>(x) > 
for x > b and <fr(x) > for x < a, which is due to the assumption that 
strict inequality holds in (16.4.6|) . We have to deal with the jump on (a, 6). 
The diagonal entries in the jump matrix on (a, b) are highly oscillatory. The 
goal of the next transformation is to turn these highly oscillatory entries into 
exponentially decaying ones. 

6.4.4 Second transformation: opening of lenses 

We open up a lens around the interval [a, b] as in Figure 16.11 and define 



/ '( _ , ) in the upper part of the lens, 

■e 2n <t> 0/ 



S = < 



1 0' 

e 2n<t> o 



in the lower part of the lens 
elsewhere. 



3.4.12) 
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Then S is denned and analytic in C \ S5 where Eg is the contour consisting 
of the real line, and the upper and lower lips of the lens, with orientation as 
indicated by the arrows in Figure 16.11 The orientation induces a +-side and a 
— side on each part of £5, where the +-side is on the left as one traverses the 
contour according to its orientation and the — side is on the right. We use S± 
to denote the limiting values of S on Eg when approaching T,g from the ±-side. 
This convention about ±-limits, depending on the orientation of the contour is 
usual in Riemann-Hilbert problems and will also be used later on. 
Then S satisfies the following RH problem: 

(1) S is analytic in C \ (R U S5). 

(2) On Eg we have the jump 5+ = S-Js where 



(3) S(z) =1 + 0(1/3) as z -¥ 00. 

If the lens is taken sufficiently small, we can guarantee that Re</> < on 
the lips of the lens. Then the jump matrices for S tend to the identity matrix 
as n — > 00 on the lips of the lens and on the unbounded intervals (—00, a) and 
(6,00). 

6.4.5 Outside parametrix 

The next step is to build an approximation to S, valid for large n, the so- 
called parametrix. The parametrix consists of two parts, an outside or global 
parametrix that will model S away from the endpoints and local parametrices 
that are good approximation to S in a neighborhood of the endpoints. 
The outside parametrix M satisfies 

(1) M is analytic in C \ [a, b], 

(2) M has the jump 




/ 1 

I e 2n<p(x) 



for x £ (a, b), 



on the lips of the lens 



(6.4.13) 



for x < a or x > b. 




for x £ (a, b). 



(3) M(z) = 1 + 0(l/,z) as z -»■ 00. 



To have uniqueness of a solution we also impose 
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Figure 6.2: Neighborhood U$(b) of b and contours and jump matrices for the 
RH problem for P 



(4) M has at most fourth-root singularities at the endpoints a and b. 

The solution to this RH problem is given by 

/ /3(z)+l3- 1 (z) f3(z)-p-\z) \ 
2 2i 
P{z)-p- l {z) P(z) + P- I (z) 



M(z) 



V 



2/ 



where 



z-b 



z — a 



7 



1/4 



z€C\[a,b], (6.4.14) 



.4.15) 



6.4.6 Local parametrix 

The local parametrix is constructed in neighborhoods Us(a) = {z \ \z — a\ < 5} 
and Ug(b) = {z \ \z — b\ < 5} of the endpoints a and b, where 5 > is small, 
but fixed. 

In Us(a) U Us(b) we want to have a 2 x 2 matrix valued function P satisfying 
the following (see also Figure HPI for the jumps in the neighborhood of b): 

(1) P is defined and analytic in (Us(a) U U$(b)) \ Eg and has a continuous 



extension to \ Us(a) U U$(b)j \ £5. 
(2) On £5 n (U s (a) U U s (b)) there is the jump 

P+=P-J S 
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Figure 6.3: Contour and jump matrices for the Airy Riemann-Hilbert prob- 
lem. 

where Jg is the jump matrix (|6. 4.13ft in the RH problem for S. 

(3) P agrees with the global parametrix M on the boundaries of Us (a) and 
Us(b) in the sense that 

P(z) = (j + O (^) \ M(z) (6.4.16) 

as n — > oo uniformly for z G dU${a) U dUg(b). 

(4) P(z) remains bounded as z — > a or z — > b. 

The solution of the RH problem for P is constructed out of the Airy function 

yo (z) = Ai(z), 

and its rotated versions 

yi(z) = uj Ai(wz), y 2 (z) = J 2 Ai(io 2 z), u = e 27ri/3 . 

These three solutions of the Airy differential equation y"(z) = zy(z) are con- 
nected by the identity 

2/0 (z) +yi(z) +V2(z) = 0. 

They are used in the solution of the following model RH problem posed on 
the contour T,a shown in Figure HT3l 

(1) A : C \ S A -> C 2x2 is analytic. 
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(2) A has jump A + = A^Ja on with jump matrix J a given by 

for arg z = 0, 



1 
-1 o, 



for arg ,z = ±2iri/3, 
for arg z = ir. 



3.4.17) 



(3) As z —7- 00, we have 



z" 1 / 4 \ 1 (\ 




-3/2- 



e 3" 



2.3/2 



2 3/2 



.4.18) 



(4) j4(z) remains bounded as z — > 0. 



The rather complicated asymptotics in (|6.4. 18[) corresponds to the asymp- 
totic formulas 



Ai(s) 
Ai'(*) 



1 



z l/4 



_2 3/2 
e 3 2 



l + (z- 3 / 2 )j , 

!--(i + 0(2 -^)), 



(6.4.19) 



as z -f with — 7r < argz < 7r, that are known for the Airy function and its 
derivative. 

The solution of the Airy Riemann-Hilbert problem is as follows 



A(z) 



2ir x < 



2/0(2) -2/2(2) 

-iy'o(z) iy' 2 (z) 

-2/1(2) -1/2(2)' 

-iyi(z) iy^(z) 

-2/2(2) 2/1(2) 

^2/2(2) - ? Vi(2) 

2/0(2) 2/1(2;) 

-^0(2) -Vl(2) 



< argz < 27r/3, 
27r/3 < arg z < ir, 
—ir < argz < — 2-7T/3, 
-27r/3 < argz < 0. 



3.4.20) 



The constants y2ir and ±i are such that det A(z) = 1 for z £ C \ E^. 

To construct the local parametrix P in the neighborhood U$(b) of b we also 
need the function 



/(*) 



(2) 



2/3 



22 



CHAPTER 6. UNIVERSALITY 




Figure 6.4: Contour and jump matrices for the RH problem for R. 



which is a conformal map from U$(b) (provided 5 is small enough) onto a neigh- 
borhood of the origin. For this it is important that the density py vanishes as a 
square root at b. We may assume that the lens around (a, b) is opened in such 
a way that the parts of the lips of the lens within U$(b) are mapped by / into 
the rays argz = ±27r/3. 

Then the local parametrix P is given in Us{b) by 

P(z) = E n (z)A(n 2 ^f(z))( en y JLA zeU s (b)\E s , (6.4.21) 



e ~™t>{*) 

where the prefactor E n (z) is given explicitly by 

Wl -i\ (nV*{z-a)V* 
nlj " v^V-^ l/\ n-V6( z _ a )-i/4 

xf (/( ^" 6))1/4 ° 1/4 V (6.4.22) 

Then E n is analytic in U$(b) and it does not change the jump conditions. It 
is needed in order to satisfy the matching condition (|6.4. 16[) on the boundary 
\z - b\ = S of U 8 (b). 

A similar construction gives the local parametrix P in the neighborhood 
Ug(a) of a. 

6.4.7 Final transformation 

In the final transformation we put 



\S(z)M(z)-\ zeC\(E s UUs(a)UU d (b)), 
\ S(z)P(z)-\ z E (U s (a) U U s (b)) \ E S . 



Then R has an analytic continuation to C\£# where is shown in Figure [67 
R solves the following RH problem. 
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(1) R is analytic in C \ 

(2) R satisfies the jump conditions R+ = R-Jr where Jr are the matrices 
given in Figure | 



(3) R{z) = 1 + Oiz- 1 ) as z — y oo. 

The jump matrices tend to the identity matrix as n — > oo. Indeed, the jump 
matrix Jr on the boundaries of the disks satisfies 

J R ( Z ) = P(z)M{zY l = I + o{ n - 1 ) 

as n — > oo, because of the matching condition (|6.4. 16|) . On the remaining parts 
of Hr we have that Jr = I + 0(e~ cn ) as n — > oo for some constant c > 0. 

Technical estimates on solutions of RH problems, see |Dei99a| . now guaran- 
tee that in this case 

R{z) = l + 0(l/n) asn^oo, (6.4.24) 
uniformly for z G C \ T,r. 

6.4.8 Proof of bulk universality 

Now we turn our attention again to the correlation kernel (|6.4.3p which, since 
n = N, we denote by K n instead of K n ^. We follow what happens with K n 
under the transformations Y h- > T t— > S. We assume that x,y G (a, b). 
From the definition ()6.4.8j) we obtain 



y/ e -nV(x)y/ e -nV(y) ,,,.,„,. . / e n(g+(a 

K n (x,y) = 2m{x _ y) (0 e^+(^/ 2 )) ^(yjT+Cz) ^ Q 

which by (non-trivial) properties of (|6.4.7p and ()6.4.9p . based on the variational 
condition (|6.4.5p . can be rewritten as 

Kn{x ' y) = 2vi{x-y) ( ° e ^ T +^ T + (*) ( J" (6A25) 

Using the transformation ([6.4. 12[) in the upper part of the lens, we see that 
(|6.4.25l) leads to 

i / p —n<f>+(x) 
Kn{x ' y) = 2m(x - y) (- en4+iy) e ~ n " +iy ^ S + ^ S +W [ e n Mx) 

(6.4.26) 

for x, y G (a, b). This is the basic formula for K n in terms of S. 
If x,y e (a + 5,b- 5), then by (|6.4.23p 

ST 1 (y)S + (x) = M~\y)R{y)- l R(x)M + {x). 
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The uniform estimate (|6.4.24p on R can then be used to show that 

S+\y)S + (x) = I + 0(x-y). (6.4.27) 

with an O-term that is uniform for x, y £ (a +5, b— S). Using (|6.4.27|) in (|6.4.26|) 
we arrive at 

ir t \ sin(m(0+(y) - <p+{x)) (*ais\ 

Kn{x,y) = -, r 1-0(1)- (6.4.28) 

vr(x - y) 

Take x* £ (a, b) fixed and let c = pv(x*) > 0. We may assume that 5 > 
has been chosen so small that x* £ (a + 5,b — 5). Replace x and y in (|6.4.28p by 
x * + cn an d x * + m respectively. Then after dividing through by cn, we have 
for fixed x and y, 

1 -r^ ( * X * y 

—K n (x* + —,x* + — 

CTl V CTl CTl 

_ sin {in {4> + (** + - *+ (x* + m +0 fl_ ] (6A29) 



7r(x — y) \n 

as n — t- oo, and the O-term is uniform for x and y in compact subsets of 
Since 4>' + (s) = iripv(s) for s 6 (a, 6), and c = /3y(x*), we have that 

in U+ fx* +y-)-<p + (x* + — )) = ir(x - y) + O ( ' ~ " 



cn' \ cn / / \ n 



and therefore the rescaled kernel (|6.4.29p does indeed tend to the sine kernel 
(I6.2,3|) as n — > oo. This proves the bulk universality. 

Remark. Letting y — > x in (|6.4.28p and dividing by n we also obtain 

-K n (x,x) = -<p' + {x) + o(- 
n tti \n 

1 



h{x)J(b-x)(x-a) + 0[ - (bv (16X91) ) 

= Pv (x) + o(^\ (bydSXS])) 

as n — )• oo, uniformly for x in compact subsets of (a, b), which proves that py 
is the limiting means eigenvalue density. 

6.4.9 Proof of edge universality 

We will not give the proof of the edge universality in detail. The proof starts 
from the representation (|6.4.26[) of the correlation kernel in terms of S. In the 
neighborhood of a and b we have by (|6.4.23[) that S = RP, where P is the local 
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parametrix that is constructed out of Airy functions, as it involves the solution 
of the Airy Riemann-Hilbert problem. 

From (|6.4.20p it can be checked that the Airy kernel (|6.2.4p is given in terms 
of the solution A(z) by 

( /- \ 

> l)A- l (y)A + (x) 



K Ai (x,y) 



1 



2iri(x — y) 



x < 



-1 l)A- l (y)A + (x) 
l)A- + \y)A + {x)i^ 
-1 l)A- l {y)A + {x) 



V 

o, 



if both x, y > 0, 
if x > and y < 0, 
if x < and y > 0, 
if both x, y < 0, 



and this is exactly what comes out of the calculations for the scaling limit of 
the eigenvalue correlation kernels K n near the edge point. 



6.5 Non-standard universality classes 

The standard universality classes (sine, Airy and Bessel) describe the local 
eigenvalue statistics around regular points. 

In the unitary ensemble ()6.3.1|) there are three types of singular eigenvalue 
behavior. They all depend on the behavior of the global eigenvalue density py . 
The three types of singular behavior are: 

• The density py vanishes at an interior point of the support. 

• The density py vanishes to higher order at an edge point of the support 
(higher than square root). 

• Equality holds in the variational inequality (|6.4.6|) at a point outside 
supply). 

In each of these cases there exists a family of limiting correlation kernels 
that arise in a double scaling limit. In (|6.3.1|) and (|6.3.3|) one lets n, N — > oo 
with n/N — > 1 at a critical rate so that for some exponent 7, the limit 

lim n 7 - l) (6.5.1) 

exists. The family of limiting kernels is parametrized by the value of the limit 
ETT) . 

Most rigorous results that have been obtained in this direction are based on 
the RH problem, and use an extension of the steepest descent analysis that was 
discussed in the previous section. Non-standard universality classes are also 
discussed in Chapters 12 and 13. 
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6.5.1 Interior singular point 

An interior singular point is a point x* where py vanishes. Varying the pa- 
rameters in V may then either lead to a gap in the support around x* , or to 
the closing of the gap. The local scaling limits at x* depend on the order of 
vanishing at x*, that is, on the positive integer k so that 

p v {x) = c(x - x*) 2k (l + o(l)) asx^-x* 

with c > 0. 

The case of quadratic vanishing (i.e., k = 1) was considered in [Ble03j for 
the critical quartic potential 

V(x) = \x" - x 2 
4 

and in [Cla06|, [Shc08j for more general V. For k = 1 one takes 7 = 2/3 in ([6.5. ip . 
The limiting kernels are parametrized by a parameter s which is proportional 
to the limit (IBToTT]) 

s = c lim ra 2//3 [ — — 1 



where the proportionality constant c > is (for the case supp(^y) = [a, b]) 

2 



(irp'yix*)) 1 / 3 y/(b - x*)(x* - a) 

The s-dependence is then governed by the Hastings-Mcleod solution q(s) of the 
Painleve II equation 

q" = sq + 2q 3 . (6.5.2) 

The limiting kernels are therefore called Painleve II kernels and we denote them 
by K pll (x, y; s), even though q(s) itself does not appear in the formulas for the 
kernels. What does appear is a solution of the Lax pair equations 

d /® 1 (x;s)\_(-4ix 2 -i(s + 2q 2 ) Axq + 2ir \($i(x;s) 
dx \$2{x; s)J ~ \ 4xq - 2ir Aix 2 + i(s + 2q 2 ) j \$ 2 (x; s) 

d_ ( $ x (x; s)\ = (-ix q\( $i(ar; s) N 
0s \$ 2 (x;s)J \q ix J \$ 2 (x; s) 

(6.5.3) 

where q = q(s) is the Hastings-Mcleod solution of (|6.5.2p and r = r(s) = q'(s). 

There is a specific solution of ([6.5.3P so that the family of limiting kernels 
takes the form 

K™(x, y; s) = + ^ s) ^ s) . (6.5.4) 

27ri(x — y) 

For k > 2, one takes 7 = 2k /(2k + 1) in (|6.5.ip . and the limiting kernels 
can be described in a similar way by a special solution of the fcth member of 
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the Painleve II hierarchy. The functions in the kernels itself are solutions of the 
associated Lax pair equations. 

The connection with Painleve II also holds for the unitary matrix model 
with a spectral singularity (I6.3.18P 

J—\detM\ 2a e - NTrV ^dM, a > -1/2. 

Z n ,N 

In the multicritical case where py vanishes quadratically at x = 0, the limiting 
kernels are associated with a special solution of the general form of the Painleve 
II equation 

q" = sq + 2q 3 — a 
with parameter a, see [Ake98l [Cla08 . 

6.5.2 Singular edge point 

A singular edge point is a point x* where py vanishes to higher order than 
square root. The local scaling limits depend again on the order of vanishing. 
If x* is a right-edge point of an interval in the support, then there is an even 
integer k such that 

p v (x) = c(x* -x) k+l,2 (l + o(l)) asx/x* (6.5.5) 

with c > 0. Here one takes 7 = (2k + 2) /(2k + 3) in (IHXTl) . The limiting 
kernels are described by Lax pair solutions associated with a special solution 
of the feth member of the Painleve I hierarchy. For k = 2 this is worked 
out in detail in [Cla07]. For general k, and assuming x* is the right-most 
point in the support, the largest eigenvalue distributions are studied in [Cla09| . 
These generalizations of the Tracy- Widom distribution are expressed in terms 
of members of a Painleve II hierarchy. 

The case k = 1 in (|6.5.5jl cannot occur in the unitary matrix model (|6.3.ip . 
although it frequently appears in the physics literature, see e.g. [D1F95], where, 
for example, it is associated with the potential 

V(x) = -—x 4 + -x 2 . 
v ' 48 2 

This potential does not satisfy (|6.3.2f) and the model (16.3. ip is not well-defined 
as a probability measure on Hermitian matrices, although it can be studied in 
a formal sense, see e.g. in Chapter 16. 

In |Fok92| and [Dui06] the polynomials that are orthogonal on certain con- 
tours in the complex plane with a weight e~ NV ( x ) t where V(x) = t^- + ^- and 
t ~ — A , are studied and the connection with special solutions of the Painleve I 
equation 

q" = 6q 2 + s 

is rigorously established. 
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6.5.3 Exterior singular point 

An exterior point x* is singular if there is equality in the variational inequality 
(|6.4,6|) , At such a point a new interval in the support may arise when perturbing 
the potential V. The limiting kernels at such a point depend on the order of 
vanishing of 

2 f ]og—*—^dnv(y) + V(x)-£ 
J \x-y\ 

at x = x* . 

In this situation the appropriate scaling is so that 

n ( n \ v 
lim — - 1 (6.5.6) 



n->oo log n \N 

exists, instead of (|6.5.ip . This special kind of scaling was discussed in |Ake97aj 



Eyn06| . Maybe surprisingly, there is no connection with Painleve equations in 



this case. In the simplest case of quadratic vanishing at x = x* the possible 
limiting kernels are finite size GUE kernels and certain interpolants, see |Ber09l 
ICla081[Mo08] . 

6.5.4 Pearcey kernels 

The Painleve II kernels ()6.5.4p are the canonical kernels that arise at the closing 
of a gap in unitary ensembles. 

Brezin and Hikami [Bre98j were the first to identify a second one parameter 
family of kernels that may arise at the closing of a gap. This is the family of 
Pearcey kernels 

K Pear (x y s) = E^KM -P'(x)q'(y) +p"(x)q(y) - sp(x)q(y) g 

x-y 

with s £ M, where p and q are solutions of the Pearcey differential equations 
p"'(x) = xp(x) — sp'(x) and q"'(y) = yq(y) + sq'(y). The kernel is also given by 
the double integral 

K J> ™(x,y;s) = I H e -^ 4 +^ 2 -™+i? 4 -!? 2 +^^ (6.5.8) 

where the contour C consists of the two rays from ±ooe i7r//4 to together with 
the two rays from to ±ooe _47r//4 . 

The Pearcey kernels appear at the closing of the gap in the Hermitian matrix 
model with external source 

JL e -nTr(V(M)-AM) dM (6 _ 5 _ g) 
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where the external source A is a given Hermitian nxn matrix. This was proved 
in |Bre98j for the case where V(x) = and A is a diagonal matrix 

A = diag(a, . . . , a, —a, . . . , —a) (6.5.10) 

n/2 times n/2 times 

depending on the parameter a > 0. For a > 1 the eigenvalues of M accumulate 
on two intervals as n — > oo and for < a < 1 on one interval. The Pearcey 
kernels (|6.5.8|) appear in a double scaling limit around the critical value a = 1, 
see also |Blc07] for an analysis of an associated 3x3 matrix valued RH problem. 
The matrix model with external source (I6.5.9P with quadratic potential has 
an interesting interpretation in terms of non-intersecting Brownian motions 
|Apt05 . 

The case of a quartic polynomial potential V(x) = ^x 4 — was analyzed 
recently in [Blellj . Here it was found that the closing of the gap can be either 
of the Pearcey type or of the Painleve II type, depending on the value of t € M. 
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